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We show that the correction-to-scaling exponents in two-dimensional percolation are bounded by 
Q < 72/91, uj = DO. < 3/2, and Ai = uui < 2, based upon Cardy's result for the critical crossing 
probability on an annulus. The upper bounds are consistent with many previous measurements 
of site percolation on square and triangular lattices, and new measurements for bond percolation 
presented here, suggesting this result is exact. A scaling form evidently applicable to site percolation 
is also found. 



In percolation, a quantity of central interest is the size 
distribution n s (p), which gives the number of clusters 
(per site) containing s sites, as a function of the site or 
bond probability p. In the scaling limit, in which s is 
large and p — p c is small such that (p — Pc)s a is constant, 
n s (p) behaves as [T] 

n s (p)~ A S -\f(B(p-p c )s°) , (1) 

where t, a, and f{z) are universal, while the metric 
constants A and B and the threshold p c are system- 
dependent. For two-dimensional systems, r = 187/91, 
a = 36/91. The constants A and B are uniquely defined 
if one assumes for example /(0) = /'(0) = 1 or /(0) = 1 
and f{z)dz = 1. 

Here we are concerned with finite-size corrections to 
Eq. Q at p c , where n s (p c ) ~ As T . It is generally hy- 
pothesized that the leading corrections to this behavior 
are of the form 

n s ( Pc )~ As- T (l + Cs- Q + ...) (2) 

where SI is the correction-to-scaling exponent. This term 
corresponds to a correction in terms of a length scale 
L ~ s 1/D of the form (1 + cL~ u + ...), where w = DSl 
and D = 91/48 is the fractal dimension. 

Studies of SI go back to the mid-1970's and are summa- 
rized in Table[l] this is an updated version of a table given 
in [9], which surveyed the work up to 1983. The values 
were found by Monte-Carlo simulation, analyses of series 
expansions, and theoretical arguments. Note that Adler, 
Moshe and Privman [5] found SI = 0.66 ± 0.07 from an 
analysis of series at p c , but found SI ~ 0.49 when studying 
the scaling behavior of the mean cluster size for p < p c , 
where corrections are proportional to (p c — p) Al and Ai 
is related to SI by A x = ((3j)tt = (91/36)0. Theoretical 
analyses of corrections to scaling were given by Aharony 
and Fisher jTUj and Herrmann and Stauffer [20] . Derrida 
and Stauffer |21j studied uj on strips, and their values 
for two different orientations extrapolate to a value of w 
somewhere in the range of 1 to 2. Bhatti et al. find Ai in 
the range w 1 — 3 for a variety of measurements related 
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to backbones of 2D percolation clusters 22 . In more 
recent work, Ziff and Bablievski [T5] found numerically 
that SI — 0.77 ± 0.02, and independently Kammerer et 
al. [19] found the identical value, SI = 0.77 ± 0.04. A 
nearby value 0.73 ± 0.02 has also been given recently by 
Tiggemann [T5] . 

In this Letter, we derive a value for SI by relating 
Cardy's recent result [23] for crossing in an annulus to 
the cluster size distribution. For the probability II(t) 
that a crossing occurs in an annulus or a cylinder, Cardy 
found 

7 7 (-l/3r)7 7 (-4/3r) = ( , ]1/2 v(3t) V (3t/4) 
[ > r ? (-l/r)r 7 (-2/3r) 1 1 ! »j(t)»/(3t/2) 

(3) 

where ry(r) q'^Un^A 1 <? n ) 

J2n=-ooh 1 ) n 1 {6n+1)2/2i is the Dedekind 77 func- 
tion, with q = e 2mT . For an annulus of outer and inner 
radii R and R\ respectively, r = (i/7r) \n{R/R{). For 
a cylinder of circumference t and length L, r = 2iL/£, 
since by conformal transformation R/Ri corresponds 
to e 27rL / l _ Note that ?j(t) can be directly evaluated in 
Mathematica using DedekindEta [r] . 

Eq. (|3| is consistent with several previous numeri- 
cal measurements of crossing. For a cylinder of as- 
pect ratio L/£ — 1/2, which corresponds to an annu- 
lus with R/Ri = e*, r = i and Eq. ^ gives n(i) w 
0.876 631451..., in agreement with the precise numer- 
ical value 0.876 657(45) found by de Oliveira, Nobrega, 
and Stauffer [SJ. For L/C = 1 or R/R x = e 27r , r = 2i 
and Eq. ^ gives n(2i) ss 0.636 454 001, which agrees 
closely with the measured values 0.636 65(8) of Hovi and 
Aharony [55], 0.638 of Acharyya and Stauffer [35], 0.64(1) 
of Ford, Hunter, and Jan [27], 0.6365(1) of Shchur [25], 
and 0.6363(3) (average) by Pruessner and Moloney [25] , 

Note that Eq. ^ implies II = 1/2 for L/t « 1.368 800, 
and the maximum in —Ti!{L/t) [the inflection point of 
n(L/f)] occurs at L/t » 0.540 652, where -II' has 
value 0.522 282. In comparison, for a system with open 
boundaries, the maximum in —Tl'(L/t) occurs at L/t « 
0.523 522 with value 0.737322 [30l E], as follows from 
Cardy's open-boundary crossing formula [32 . 
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TABLE I. History of determinations of fi, u — DO. = (91/48)fi, and Ai = Q/a = (91/36)f2. Numbers in parentheses represent 
errors in last digit (s), and are shown on original values only. 



year 


author 


method 




w 


Ai 


1976 


Gaunt & Sykes [2] 


series 


0.75(5) 


1.42 


1.90 


1978 


Houghton, Reeve & Wallace [3] 


field theory 


0.54-0.68 


0.989-1.28 


1.32-1.71 


1979 


Hoshen et al. p[] 


MC 


0.67(10) 


1.27 


1.69 


1980 


Pearson [5] 


conjecture 


64/91 « 0.703 


1.333 


1.778 


1980 


Nakanishi & Stanley [B] 


MC 


0.6 < Q < 1 






1982 


Nienhuis [7] 


field theory 


96/91 w 1.055 


2 


2.667 


1982 


Adler, Moshe & Privman [8] 


series 


0.5 


0.95 


1.26 


1983 


Adler, Moshe & Privman [9] 


series 


0.66(7) 


1.25 


1.67 




53 


series p < p c 


0.49 


0.93 


1.24 


1983 


Aharony & Fisher [TU1 [TT] 


RG theory 


55/91 w 0.604 


55/48 « 1.15 


55/36 « 1.53 


1984 


Margolina et al. [Ill 112] 


MC 


0.64(8) 


1.21 


1.62 






series 


0.8(1) 


1.52 


2.02 


1985 


Adler [13] 


series 


63(5) 


1.19 


1.59 


1986 


Rapaport [14] 


MC 


0.71-0.74 






1998 


MacLeod & Jan [T3] 


MC 


0.65(5) 


1.23 


1.64 


1999 


Ziff and Bablievski [16] 


MC 


0.77(2) 


1.46 


1.95 


2001 


Tiggemann [T7] 


MC 


0.70(2) 


1.33 


1.77 


2007 


Tiggemann [18] 


MC 


0.73(2) 


1.38 


1.85 


2008 


Kammerer, Honing & Franosch [19] 


MC 


0.77(4) 


1.46 


1.95 


2010 


this work 


theory 


72/91 « 0.791 


1.5 
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Expanding the rj functions in ([3| for large L/£, we find, 

U(q) = yf<f /48 (l ~ <f /2 + f - q V2 + 2q 4 - if' 2 -..) 

(4) 

where q = R x /R = e -^L/t ag in j^gj Thugj for thc 
annulus, we have 

+ (4r-(ir-(ir--}» 

and for the cylinder, we have 

n ( L/*) = y| e -(5/24)^|: _ e - 3 ^/^ 

+ e -^i/< _ e -™7* + 2e - 8 ^/^ - ... | (6) 

Shchur |28j has verified the leading term numerically, 
finding 0.65448(5) for the exponent 5tt/24 ps 0.654498, 
and the intercept of the asymptotic line in his Fig. 
6 is consistent with the predicted coefficient \/3/2 w 
1.224 745. Preussner and Moloney [33] also measure this 
coefficient and find 1.2217(4) for bond percolation and 
1.2222(4) for site percolation on a square lattice, where 
the errors bars represent statistical errors of a single set 
of systems. 

For the annular result, we can imagine that the system 
is actually infinite with an inner boundary of radius Ri, 
then, H(R/Ri) is the probability that a cluster connected 



to that inner boundary has a maximum radius greater 
than or equal to R. That is, if, in the infinite system, the 
cluster connected to the center extends beyond a circle 
of radius R, then in the annulus there will be a crossing 
cluster between the two circles of radii Ri and R, and 
these two events will occur with the same probability. 
Thus II gives a measure of the size distribution of the 
clusters connected to the inner circle, where the size is 
characterized by the maximum cluster radius. Here, we 
relate this to thc cluster size distribution by associating 
the inner radius to the discreteness of the lattice. 

Given a size distribution n s , the probability that an 
occupied site is connected to a cluster of size greater than 
or equal to s is given by: 

°° poo 

p >s(p) - s ' ns ' ~ / s ' ns,ds ' ( ? ) 

Using |2]) for n s (p c ), we thus find that, at p c , 

P> s (p c )~A' s 2 - T (l + B> s - n + ...) (8) 

where 2 — r = — j3a = —5/91 in 2-d. Because critical 
percolation clusters are fractal, s are R are related by 

8 ~ so(R/e) D (9) 

where D is the fractal dimension, sq is a constant, and e 
is of the order of the lattice spacing and represents the 
lower size cutoff of the system, similar to a boundaryt 
extrapolation length [33]. Putting ([9| into (6), assuming 
e = Ri , we find that the probability a cluster has a radius 
greater than R is given by 

P> R ( Pc ) ~ aR D ^(l + bR- nD + ...) (10) 
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where a and b are constants. By hyperscaling D(2 — r) = 
D — d = —fi/v = —5/48. Comparing this with (j5j), we 
see that 17£> = w = —3/2, or 



17 = 3/ (2D) = 72/91 



(11) 



implying also Ai = uuj = 2. 

Alternately, if we put (R/Ri) = (s/s ) 1/£l into (JHJ), we 
find 



P> s (Pc) = 



3 / s 



2 Vso 



1 - 



so 



•so 



.(12) 



which gives the higher-order corrections also. 

In deriving ( 12 ) we have ignored finite-size corrections 
to ([9]). Say we have to next order 



s~s (R/e) D (l + cR- x + ...) 



(13) 



then this will lead to a term of order (s/sq)~ x ' d in 
the expansion of P> s (p c ), and where x/D lies within 
the exponents 72/91, 96/91, ... will determine its im- 
portance. It is of course possible that x/D < 72/91, 
in which case that would be the dominant correction. 
There have been studies in the past on the finite-size 
corrections to the radius of gyration [34], but not to our 
knowledge to the maximum radius with respect to an 
arbitrary point within a cluster, which is needed here. 
Another source of finite-size corrections are discussed by 
Fisher and Aharony (see footnote 18 of [H]) which imply 
a correction 17 = 55/91 = 0.604. 

Our result for 17 is consistent with the very first deter- 
mination 17 = 0.75 ± 0.05 in 1975 by Gaunt and Sykes, 
and many of the subsequent results including the latest 
numerical values. Some previous results gave lower val- 
ues, such as the series results of [5J, which were based 
upon studying the scaling of quantities like the mean 
cluster size away from p c . Perhaps these lower values 
were due to analytic corrections which are relevant away 
from p c . 



In Eq. ( 12 1, we see that the next-order correction term 
is s~ n2 with 17 2 = 96/91. This interestingly is exactly the 
value of 17 proposed by Nienhuis [7] . The closeness of this 
exponent to 17 would make its determination numerically 
difficult, and in any case higher-order corrections in the 
relation ^ between s and R might mask this term. 

Assuming the other correction terms are small, Eq. 



(12 1 implies that P> s should be a universal function of 
s/sq, where sq varies from system to system. This may 
explain an observation made in |16| that the quantity 
17(est) defined by 



fi(est) = — log 2 



C s — C. 



a/2 



c 



a/2 



c 



s/4 



(14) 



where C s = s T 2 P> S , when plotted vs. In s, appears to be 
a universal curve for site percolation on the square and 




FIG. 1. Plot of — S7(est) as a function of ln(s/so) for s = 2, 
4, ... 1024. The two upper curves are a replotting of the 
data given in Fig. 6 of Ref. J16 , using so = 0.25 (square 
lattice) and so = 0.13 (triangular lattice) to line the data 
up with the theoretical prediction — the first three terms of 
(12 1. No vertical rescaling was done. (Note: there was a 
minus sign missing in the label of the vertical axis in Ref. 
[16] .) The two lower curves are bond percolation on the square 
and triangular lattices, using so = 0.25 and 0.5, respectively. 
[Note: data will be updated for published version] 



triangular lattice, except for a horizontal shift. 17(est) 
equals 17 if the correction to scaling has only one term as 
in {[sj) , and otherwise can be used to estimate 17 by taking 
s relatively large. The shift that was needed in [16] to line 
up the data just reflects the difference in s between the 
two lattices, because of the logarithmic scale in the plot. 
In Fig. [T] we have replotted the data of Ref. 16] along 
with the results of using the Q(est) that follow from ( 12 ) 



(using the first three terms), and adjusted sq to match 
the theoretical prediction. The behavior can be seen to 
match the theory very well, with only a single adjustable 
parameter s - In Ref. [TB] , we generated s» 6- 10 9 clusters 
up to size s = 1024 for each system, to obtain these data. 

Here we have generated additional data for bond per- 
colation on square and triangular lattices, in which we 
characterize s by the number of sites wetted by the clus- 
ters. (Evidently, these systems have not been studied in 
the past to find 17.) We generated ~ 10 11 clusters with 
s < 1024 for each lattice, using the R(9689) random- 
number generator of |35j . The corresponding 17(est) are 
also plotted in Fig. [I] For both lattices 17 (est) approaches 
~ 0.8 from above for larger s, implying that (11 1 is again 
the leading correction exponent, but for smaller s the 



data do not follow the behavior implied by (12 1, imply- 



ing that for bond percolation, there are other significant 
finite-size corrections. While ( 12 1 implies deviations from 



the leading behavior are negative, for bond percolation in 
fact those deviations are positive. Similar positive finite- 
size effects are seen in the closely related SIR epidemic 
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model [3"fJ] . 

It appears that fi(est) gives a lower bound to f2 with 
site percolation, and an upper bound with bond percola- 
tion, both bracketing 0.79. Evidently larger clusters are 
needed for fi(est) to get close to that value, but statis- 
tical errors limit the cluster size that can practically be 
used. (The difficulty in reaching this asymptotic regime 
for the corrections might explain the generally low val- 
ues of found in the past.) We have also carried out 
an analysis of the exact n s (p c ) for site percolation on the 
square lattice, based upon the known series ([37JI3E]) for 
n s (p) for s < 22. Results are lead to similar conclusions 
and will be discussed elsewhere. 

We can find the behavior of n s from ( 12 ) as follows: 



n s {pc) = — (P>s - P>(s+i)) = 



5p c 
91s 2 



1 - 



77 



48 
91 



3 / * 
2 \s 
101 / s 
5 \s 



(15) 



where the factor of Vc is included for site percolation only 



PQ. Thus, for Eq. (lj), we have A = (5p c /91) 



5/91 



and B = (77/5)sq 2 ^ 91 , and here we have picked up the 
analytic term s~ x . Likewise, 



n>s(p c 



5p c 
96s 



so 



44 



so 



(16) 

For example for site percolation on the triangular lattice 
(p c = 0.5), we found that sq — 0.13 (see Fig. [T]), which 
implies a coefficient of 0.0285 to n> s (p c ). This value 
agrees closely with the value « 0.028 found in [TT| (their 
Fig. 2.) 

In conclusion, we found that the behavior of crossing 
on an annulus implies the correction-to-scaling exponent 
(111 to the size distribution n s , which seems to represent 



the dominant term for the systems we have studied, in 
comparison to other possible finite-size terms that might 
-appear. For future work, it would be interesting to study 
ts(Pc) on different lattices, and to explore directly the re- 
lation between R and s andfind higher-order corrections 
to @. 

The author thanks Amnon Aharony and Joan Adler 
for providing comments and references, and especially 
thanks Dietrich Stauffer for many helpful suggestions. 
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